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A  Ray  Trace  Program  for  Line-of-Sight 
Microwave  Communication  Links 


1.  INTRODUCTION 

An  ctT'ift  was  initiaifd  in  S<'(>t<’nib('r  1981  to  obtain  high  resolution  vertical  pro¬ 
files  of  the  nderowavi'  index  of  tad'raetion  eonlemporaneousl/  with  link  dela/ 
nu'asurenK.nts  on  line-of-sight  (t.OS)  eomni unication  links.  The  vertical  profiles  of 
the  index  of  ((Traction  were  obtained  using  a  portable  balloon  borne  radiosonde 
systoni.  Those  systems  are  capable  of  measuring  the  refractive  index  of  about  every 
a  sec  which  in  our  application  means  about  10  m  in  altitude.  These  soundings  could 
be  taken  as  close  in  time  as  iO  t  'in  and  at  more  than  one  location  along  the  link. 

An  atmospheric  propagation  model  was  sought  that  would  take  into  account 
I  (implex  terrain  and  horizontal  gradients  in  refractivity.  These  constraints  together 
with  the  high  freriuem  ies  (10-15  GHz)  made  a  ray  trace  model  seem  the  most  attrac¬ 
tive  alternative.  Although  we  used  an  existing  ray  trace  model  for  some  early  in¬ 
sights,  it  was  evident  we  needed  a  model  tailored  to  our  specific  needs.  This 
report  discusses  the  capabilities  and  mechanics  of  the  model  developed. 


(Heceived  for  publication  3  October  1984) 

1.  Wasky,  R.P.  (1977)  A  Geometric  Optics  Model  for  Calculating  the  Field 

Strength  of  Klectromagnetic  Waves  in  the  Presence  of  a  Tropospheric  Duct, 
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2.  RAY  TRACE  MODEL  DESCRIPTION 


The  computer  program  developed  was  written  in  Fortran/ V  for  use  on  the 
AFGL  Cyber  750  computer.  Because  it  is  modular,  it  is  easily  modified  for 
application  to  other  ray  trace  problems.  It  is  written  specifically  for  line -of-sight 
communication  links  and  uses  meteorological  data  from  radiosonde  measurements. 

One  version  accepts  meteorological  soundings  at  several  locations.  It  re¬ 
quires  information  on  the  vertical  profile  of  refractive  index,  the  height  of  the 
antennas,  and  the  earth  contour  directly  under  the  transmission  beam.  The 
refractive  data  are  provided  at  "significant"  levels  determined  in  the  same  general 
manner  that  standard  radiosonde  data  "significant”  levels  are  selected.  The  earth 
contour  data  are  presented  as  a  sequence  of  points  (2 -dimensional)  where  the  earth's 
surface  slope  changes  significantly.  These  points  are  nodes.  The  program  uses 
ihe  refractive  index  data  to  calculate  ray  trajectories  and  the  earth  contour  to 
determine  shielding  effects  or  reflections. 

The  model  has  two  basic  versions.  One  is  designed  to  have  a  single  vertical 
profile  of  refractive  index  that  assumes  the  atmosphere  is  horizontally  homogeneous. 
The  other  can  ac'-ept  several  profiles  along  the  link  and  then  interpolates  refractive 
structure  between  profiles.  If  the  antennas  are  outside  the  last  radiosonde,  the 
model  assumes  horizontal  homogeneity  in  the  end  portions.  There  are  several 
(  ontrol  parameters,  including  the  maximum  and  minimum  initial  ray  angle,  the 
angular  spacing  of  the  rays,  and  whether  the  operator  wants  no  earth  reflection, 
wants  specular  earth  reflection,  or  wants  all  reflected  rays  directed  to  the  receiving 
antenna. 

The  calculated  traje  tories  are  plotted  on  a  flat  earth  profile  with  a  vertical 
scale  enlarged  58  times  relative  to  the  horizontal  scale.  The  time  of  arrival  of 
I  [lose  ravs  that  reach  the  end  of  path  are  presented  in  a  second  graph  as  a  function 
if  altituile.  The  time  of  arrival  graph  shows  the  times  relative  to  the  time  of 
arrival  of  the  rav  with  the  shortest  propagation  time.  A  third  graph  gives  the  angle 
•if  arrival  of  these  ravs. 

Figures  1  through  4  are  presented  to  demonstrate  the  various  options  and 
\ersions  if  the  model  available  to  the  operator.  Data  from  a  communication  link 
in  Cermanv  over  the  Rhine  River  on  4  Nov  I'tSl  were  used  for  this  example. 

Figure  la  shows  the  version  of  the  program  designed  for  one  refractive  index 
input  profile  .  This  refractive  structure  profile  is  shown  by  the  pair  of  lines  at  a 
distam  e  of  about  67  km  along  the  path.  The  line  that  is  exactly  vertical  indicates 
the  location  on  the  link  where  the  radiosonde  was  released.  The  second  line  is  the 
refractive  inrlex  data  profile  (in  m  units)  derived  from  that  radiosonde  (see 
Appendix  A  for  ni  unit  definition).  The  m  scale  is  shown  at  the  top  of  the  plot.  The 
altitude  is  given  in  meters,  the  distance  in  km.  All  rays  emanate  from  the  antenna 
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on  the  left  hand  end  of  the  link.  Only  the  maximum  and  minimum  rays  are  plotted 
near  the  antenna  because  the  density  of  the  rays  would  cause  the  plotting  paper  to 
saturate  and  tear.  This  example  is  for  the  case  where  no  reflection  has  been 
selected  so  all  rays  that  strike  the  ground  are  truncated.  The  ground  is  represented 
by  the  height  profile  at  the  bottom  of  the  plot. 

Figure  lb  shows  the  relative  arrival  times  (delays)  of  those  rays  that  arrive 
at  the  end  of  the  link  and  have  neither  struck  the  ground  nor  gone  above  the  maxi¬ 
mum  plotting  height.  Figure  Ic  shows  the  angle  of  arrival  of  the  same  set  of  rays 
as  those  in  Figure  lb.  Both  the  arrival  time  and  angle  are  shown  at  the  altitude  of 
the  ray  at  end  end  of  the  link.  If  a  horizontal  line  is  drawn  from  lb  and  Ic  to  la 
the  ray  corresponding  to  the  plotted  arrival  time  and  delay  can  be  identified.  This 
is  done  with  a  dashed  line  at  a  height  of  800  m  for  illustrative  purposes. 

Figure  2  shows  the  same  case  as  shown  in  Figure  1,  using  the  option  that  allows 
specular  reflection.  To  keep  the  plot  (Figure  2a)  from  getting  too  confusing, 
plotting  of  t:  e  following  rays  have  been  suppressed:  reflected  rays  that  go  out  of 
the  plotted  area  at  the  top,  rays  that  have  a  reflection  angle  of  greater  than  5’,  and 
rays  that  re -intersect  the  earth  (that  is,  a  second  reflection)  are  suppressed  at 
the  original  intersection.  In  the  time  delay  plot  (Figure  2b)  and  the  angle  plot 
(Figure  2c)  the  reflected  rays  that  reach  the  end  of  the  link  are  indicated  by  the 
he.xagonal  symbols. 

Figure  ,3  shows  the  same  case  using  the  option  that  directs  all  reflected  rays 
toward  the  receiving  antenna.  This  option  adds  additional  rays  over  the  number 
determined  by  the  maximum  and  minimum  angle  and  angular  increment.  These 
additional  rays  are  chosen  by  directing  rays  toward  the  inflection  points  used  to 
describe  the  earth's  contour,  which  are  then  reflected  toward  the  receiver.  Only 
those  ravs  that  are  unshielded  in  both  legs  of  the  journey  are  plotted.  This  type  of 
reflected  ray  is  depicted  with  a  triangle  symbol  in  the  delay  (Figure  3b)  and  angle 
(  Figure  3c )  plots. 

Figure  4  illustrates  the  version  of  the  profile  designed  for  input  from  more 
than  one  radiosonde.  In  this  case  two  soundings  have  been  used,  one  at  42  km  and 
me  at  67  km.  The  rav  traiectorv  is  computed  using  only  the  42  km  sounding 
(horizontal  hrimrjgemdt v )  up  to  42  km,  refractive  data  between  42  and  67  km  are 
I  ompiled  bv  linear  interpolation,  and  beyond  67  km  only  the  sounding  at  67  km  is 
used.  An  earlier  version  using  extrapolation  in  the  outside  zones  was  tried  but 
resulted  in  .erv  unrealistic  refra<'tive  gradients  and,  consequently,  ray  traces. 

This  version  also  has  the  reflec  ti<jn  jption.s  shown  for  the  single  profile  case. 


3.  PROGRAM  DESCRIPTION 


A  general  flow  diagram  of  the  program  is  shown  in  Figures  5a,  5b,  5c,  and  5d. 
Figure  5a  shows  the  main  part  of  the  program  RAYSUB.  Figure  5b  shows  the 
subroutine  RAYTRACE.  Figure  5c  shows  the  subroutine  SEARCH.  Figure  5d 
shows  the  subroutine  RAYINC. 

The  main  program,  RAYSUB,  does  most  of  the  housekeeping  chores  and  has  two 
principal  parts.  The  first  takes  care  of  the  evenly  spaced  fan  of  rays  emanating 
from  the  transmitter.  The  second,  when  selected,  directs  rays  at  each  of  the  nodes 
used  to  describe  the  earth  profile  and  then  redirects  them  to  the  receiver.  The 
first  task  is  accomplished  bv  computing  the  initial  angle  and  elevation  for  each  ray 
and  calling  the  RAYTRACE  subroutine.  The  second  task  is  accomplished  by  using 
the  SEARCH  routine  to  calculate  the  initial  angle  to  hit  the  earth  at  a  node,  then 
R.AYTRACE  to  get  there,  then  SEARCH  again  for  the  initial  angle  t(j  hit  the  rec  eiver 
and  finally  RAYTRACE  again  to  get  there.  Other  functions  of  RAYSUB  are  house¬ 
keeping  tasks.  It  reads  the  input  refractive  data,  converts  it  to  m  units  and 
extrapolates  it  down  and  up  to  make  sure  refractivity  is  defined  in  the  complete 
space.  It  deletes  reflected  rays  that  are  shielded  and  provides  the  plotting  routines. 

The  R.AYTRACE  subroutine  (Figure  3b)  oversees  each  individual  raytrace.  It 
requires  the  starting  increment  index,  the  type  reflection  desired  (KBNCE),  the  ray 
index,  (1).  whether  to  ignore  the  earth  (ISVVICH),  and  the  position  along  the  link  to 
terminate  the  ray  (.\POS2), 

R.AYTR.ACE  selects  the  incremental  distance  to  be  used  in  computing  the  tra- 
jectorv.  It  uses  the  lesser  of  one  two-hundreth  of  the  complete  path  length  or  the 
distance  to  the  next  surface  node  (in  the  multiple  sounding  version  it  also  considers 
the  distance  to  the  next  sounding).  It  uses  RAYINC  to  do  each  increment  of  the 
trajectory  and  provides  RAYINC  with  starting  angle,  height,  increment,  and  whether 
tc  ignore  the  earth  surface.  If  the  return  from  RAYINC  shows  the  ray  to  intersect 
the  earth.  RAYTRACE  terminates  the  ray  (KBNCE  0),  or  calculates  the  angle  for 
specular  reflection  (KB.NCE  1)  and  continues,  or  calculates  the  angle  needed  to 
reflect  to  the  antenna  (KBNCE  2)  using  SEARCH  and  continues. 

The  SEARCH  (Figure  5c)  subroutine  is  an  iterative  program  for  aiming  rays 
from  one  point  to  another.  It  estimates  the  initial  angle  from  the  geometry,  uses 
RAYINC  to  calculate  the  path,  then  corrects  the  initial  angle  estimate  using  the 
error  (distance  ray  missed  antenna).  It  ignores  intersection  with  the  earth 
(ISVVICH  -1)  since  it  is  only  used  to  calculate  an  initial  angle.  Since  not  all  paths 
are  possible,  (for  example,  there  may  be  a  radio  hole),  it  aborts  after  nine  attempts 
if  the  ray  is  more  than  3  m  from  the  intended  target. 


The  RAYING  (Figure  5d)  subroutine  computes  the  incremental  ray  trajectories 
(Appendix  A).  This  routine  receives  a  height  and  an  angle  at  the  beginning  of  the 
increment  and  calculates  the  height  and  angle  at  the  end  of  the  increment.  If  this 
second  height  indicates  it  has  passed  a  '•efractively  significant  altitude  it  deter¬ 
mines  the  height  and  distance  of  the  intersection  of  the  ray  with  that  altitude. 

After  this,  if  selected,  it  tests  to  see  if  it  would  intersect  the  earth's  surface.  If 
it  would  it  finds  the  point  of  intersection.  RAYING  returns  the  height,  angle  and 
horizontal  distance  to  the  nearest  intersection,  if  any,  inside  the  original  increment. 
It  also  informs  the  RAYTRAGE  of  the  type  of  termination  (KTYPE;  1-normal; 

2 -intersect  significant  altitude;  3 -intersect  earth). 

The  ray  trace  model  described  is  being  used  together  with  fine  scale 
meteorological  profile  data  to  analyze  performance  on  active  LOS  communication 
links.  These  tests  were  performed  in  conjunction  with  channel  probe  measure¬ 
ments  of  time  of  arrival,  and  in  some  cases,  angle  of  arrival  measurements.  The 
model  generates  the  time  of  arrival  and  angle  of  arrival  plots,  which  allow  direct 
comparison  to  both  of  the  direct  channel  measurements.  Because  of  the  modular 
structure  of  the  computer  program  used,  it  is  easy  to  adapt  it  to  do  special  tasks. 

It  has,  for  example,  been  used  to  see  the  effect  of  coupling  into  a  duct  that  is  too 
high  for  "normal"  ducting  and  is  being  adapted  to  analyze  troposcatter  radio 
propagation. 
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Figure  5c.  SEARCH  Subroutine 


Figure  5d.  RAYING  Subroutine 


Appendix  A 

Computation  of  Ray  Trajectories 


The  angle  that  a  straight  line  (ray)  makes  with  any  line  normal  to  the  local 
vertical  in  a  circular  coordinate  system  (Figure  Al)  can  be  described  by 
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The  law  of  refraction  for  a  nonhomogeneous  dielectric 
be  curved  according  to 
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states  that  the  rays  will 
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Al.  Kern,  D.  E.  (1951'  Propagation  of  Short  Radio  Waves,  McGraw  Hill,  New  York. 


where  R  is  the  radius  of  curvature,  n  is  the  refractivity  of  the  dielectric  (air)  and 
the  gradient  of  refractivity  normal  to  the  ray  trajectory.  The  result  of  this 
curvature  on  6  is  analogous  to  the  effect  of  the  curvature  of  the  earth,  with  the 
opposite  sense,  that  is,  a  positive  R  results  in  a  reduction  in  the  angle.  Hence,  if 
we  use  the  same  small  angle  approximation,  |  r|  is  of  the  same  order  as  r  or 
greater,  we  can  write 
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or  using  Eq.  (A2) 
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where  M  is  the  ijiodified  index  of  refraction  given  by 
M  =  i\  +  157  z  . 


An  analytical  expression  for  the  height  Zj  of  the  ray  above  the  surface  of  the  earth 
can  be  obtained  from  the  MacLaurin  series  expansion  and  Eq.  {A3) 
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but 


Z'  =  tan  6  ~  9 
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Equations  (A3)  and  (A5)  are  sufficient  to  calculate  the  ray's  trajectory  given  the 
initial  angle  and  the  gradient  of  the  refractive  index.  It  should  be  noted  that  small 
angle  approximations  were  used  in  several  places.  The  intended  use  of  this  tech¬ 
nique  is  for  line-of- sight  communication  links  where  the  maximum  angles  are  about 
0.  3  deg.  The  product  using  this  technique  has  been  checked  against  other  ray 
trace  techniques  that  do  not  make  these  approximations  and  no  significant  differences 

have  been  noted  in  angle,  height,  or  relative  propagation  times.  Similar  treatment 

A2  A3 

can  be  found  in  Millington  and  a  more  rigorous  development  in  Dougherty. 


A2.  Millington,  G.  (1957)  The  concept  of  constant  radius  of  the  earth  in  tropospheric 
propagation,  Marconi  Rev. ,  _2p(No.  l):79-93. 

A3,  Dougherty.  H.  T.  (1981)  Electromagnetic  Wave  Trajectories  at  VHF  and  Higher 
Frequencies,  NTIA-CR-81-14. 
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